Crossover behavior of the J\-Ji model in a staggered magnetic field 
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The ground states of the S = | Heisenberg chain with the nearest-neighbor and the next-nearest- 
neighbor antiferromagnetic couplings are numerically investigated in a staggered magnetic field. 
While the staggered magnetic field may induce the Neel-type excitation gap, and it is characterized 
by the Gaussian fixed point in the spin-fluid region, the crossover to the behavior controlled by the 
Ising fixed point is expected to be observed for the spontaneously dimerized state at finite field. 
Treating a low-lying excitation gap by the phenomenological renormalization-group method, we 
numerically determine the massless flow connecting the Gaussian and Ising fixed points. Further, 
to check the criticalities, we perform the finite-size-scaling analysis of the excitation gap. 
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Theoretically and experimentally, there have been in- 
tensive investigations on critical phenomena observed in 
low-dimensional quantum spin systems. In particular, re- 
cent investigations on this subject have focused not only 
on critical fixed points and their neighboring regions, but 
also on the global behaviors of the renormalization-group 
(RG) flows connecting them. For one-dimensional (ID) 
quantum systems (also for 2D classical systems) , besides 
the exact solutions available for some cases, the confor- 
mal field theory (CFT) provides the most efficient way to 
characterize the fixed points, where the values of the cen- 
tral charge c specify their universality classes. Further, 
with respect to the RG flows observed in the continuum 
and unitary models, Zamolodchikov's c-theorem serves 
for the explanations of their general properties .u Since in 
investigations of the quantum spin chains and interact- 
ing electron systems, the c = 1 CFT has close relevance 
to their criticalityp it is very important to understand 
its instability and a crossover to behaviors controlled by 
other critical fixed points. 

In this paper, we shall study the ground states of 
the S = I Heisenberg chain with the nearest-neighbor 
and the next-nearest-neighbor antiferromagnetic cou- 
plings (the so-called J\-J 2 model) in a staggered mag- 
netic field; the Hamiltonian to be considered is given by 
H = H x +H 2 with 



Hx = J2 (2JiSj ■ S j+1 + 2J 2 S 3 ■ S j+2 ) , (1) 

j'=i 

L 

H 2 =j2-u-iys^, (2) 

3=1 

where SJ is the vth component of the spin operator on 
the jth site Sj and couplings J\, J 2 > (we take J\ as 
the energy unit in the following) . We assume the periodic 
boundary condition S^ + i = Si and an even number of 
L. 

In the zero-field case (h s = 0), there are two special 
points where the exact ground states have been known: 
at J 2 — where the Bethe ansatz solution is available, the 



system is in the spin-fluid phase described by the Gaus- 
sian fixed point (c = 1 CFT); at J 2 — 0.5 (the Majumdar- 
Ghosh point) the direct products of spin singlet pairs 
formed either on bonds (2k, 2k + 1) or (2k + 1,2k + 2) 
become twofold degenerated -ground states (i.e., sponta- 
neously dimerized states), UuU and a finite excitation gap 
exists there.tl According to HaldaneQ, and Kuboki and 
FukuyamaJ3 the gap formation caused by the frustration 
can be well described by the quantum sine-XJordon model 
obtained via the bosonization procedure;^ the effective 
Hamiltonian of Eq. ([!]) is given as 
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where the bosonic operator 9 is the dual field of (j> satisfy- 
ing the commutation relation d y [4>(x), 9(y)\ — md(x-y). 
K and v are the Gaussian coupling and the spin- wave ve- 
locity, andj^w, stands for the spin Umklapp scattering bare 
amplitudellS Although the property of this effective con- 
tinuum model has been well understood and actually the 
<70 term may become relevantjii for the determination of 
the fluid-dimer transition point J 2 , numerical treatments 
of Eq. (0) were required. By carefully investigating the 
effect of the marginal term on the critical fixed point 
through lower-energy excitation levels observed in finite- 
size systemSjXlkamoto and Nomura precisely determined 
J 2 ~ 0.2411.E3 While it is known that the excitation gap 
also exists in the case of J 2 > 0.5 Jl3 we, in the following 
discussion, restrict ourselves to the region < J 2 < 0.5 
for simplicity. 

In the case of nonzero field [this is approximately re- 
alized in the quasi-lD antiferromagnets with the_alter- 
nating |gyromagnetic tensors, e.g., Cu-benzoateO and 



Yb4As3li3 and the field may be also generated as the p.- 
trinsic one originated from the Neel ordered sublatticelij] 
since the bosonized form of Eq. (0) is given as 
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and the scaling dimension of this perturbation teem is 
Xi = | = 2 — x 2 = 3/2) on the fixed point,EZl the 

second-order phase transition occurs for systems in the 
spin-fluid region Ji < Jj* with the divergence of correla- 
tion length £ oc h s 2 ^ 3 (aside from the logarithmic correc- 
tion due to the g<j, term). On the other hand, for J2 > J| , 
due to the relevant g^ term, the crossover of the transi- 
tion driven by the staggered magnetic field occurs. Since 
the dimer gap may survive in a weak field, the critical 
point h*(J2) takes nonzero values depending on J 2 , and 
further the.-universality class is changed. Recently, Fab- 
rizio et aL,Ea on the basis of the double- frequency sine-, 
Gordon (DSG) theory given by Delfino and Mussardo,lL3 
argued that the system on hg(J 2 ) is renormalized to the 
Ising fixed point with c = ^ in accord with the "down- 
hill" condition of the c-theorem, so the criticality in the 
vicinity of this line is related to the divergent correlation 
length of the form £ cx [h s — h*(J2)]~ 1 (explained below) 
described by the </? 4 theory. 

From the viewpoint that the Gaussian fixed point is 
perturbed by two relevant operators of g$ and h s , we 
can qualitatively estimate /i*(J2) around the point ac- 
cording to the crossover argument.E3 However, to evalu- 
ate its precise value especially near the Ising fixed point, 
a numerical treatment of H should be required. Here, 
it should be noted that the criticality on the Neel-phasc 
boundaries in the S = 1 bond-alternating XXZ chain also 
belongs to the Ising universality and they were precisely 
determined by the numerical method.EIJ Thus, we shall 
perform our calculations by analyzing the lower-energy 
excitations observed in the finite-size systems. Let us 
first focus our attention to the excitations in the, s j pip - 
fluid region described by the Gaussian fixed point ,E3lH3 

O x = V2sin\/20, (5) 

2 = V2cos\/2>, (6) 

3 = exp(±h/20). (7) 

0i and 2 denote dimer and Neel excitations, respec- 
tively, while 03 is the doublet excitation changing an 
amount of the total spin. According to the finite-size- 
scaling (FSS) argument based on CFT, corresponding- 
energy levels AEi for these operators (taking the ground- 
state energy as zero) are expressed by the use of their 
scaling dimensions AEi 2-KvXi/LH For the 

case of h s — 0, we can calculate AEi according to the 
level-spectroscopy methodJ£3 where discrete symmetries 
of the lattice Hamiltonian (see below) are utilized to spec- 
ify excitation levels. On the other hand, for the case 
of h s ^ 0, the usable symmetry becomes lower, and 
more importantly, the universality class is changed to 
the Ising one on the line h*(J2) so that we should em- 
ploy other criterion to characterize the levels. Here, we 
will use the so-called-XQZ-IR (ultraviolet-infrared) oper- 
ator correspondence Ji§E3 i.e., along the critical RG flow 
the operators on the Gaussian fixed point (UV) are trans- 
muted to those on the Ising fixed point (IR) as 




FIG. 1. The L and hs dependences of LAE( J 2 , /i s , L). 
From left to right, the next-nearest-neighbor coupling 
J 2 =0.32, 0.40, 0.46, and 0.5. Data for systems of L = 18 
(crosses), 20 (triangles), 22 (squares), 24 (diamonds), 26 (cir- 
cles), and 28 (double circles) are plotted with the fitting 
curves. 

d-M. 2 ->I+e, (8) 

where fi is the disorder field (Z2 odd) and e is the en- 
ergy density operator (Z2 even) with scaling dimensions 
Xfj, = \ and x e = 1, respectively. According to this corre- 
spondence, we can obtain a relevant excitation in nonzero 
field by taking the limit h s \ and assigning it to one 
of well-characterized states. From Eqs. (^) and (@), the 
staggered magnetic field plays a role of the "thermal" 
scaling variable on hg(J 2 ), and thus the critical exponent 
1/v = 2 — x e = 1 as already mentioned (see Ref. [18]) . On 
one hand, since the "magnetic" excitation /i stemming 
from the dimer excitation provides lower energy (i.e., the 
most divergent fluctuation), we will thus focus attention 
on it for the determination of h*(J2). 

Now, we shall explain the numerical calculation proce- 
dure. Since the system is massive unless it is located on 
the critical line, the phenomenological renormalization- 
group (REG) method is expected to work efficiently for 
our aimj£3 We numerically solve the PRG equation for 
the systems of L and L + 2, 

(L + 2)AE(J 2 , hs, L + 2) = LAE{J 2 , h s , L) (9) 

with respect to h s for a given value of J2. Since the 
equation can be satisfied by the gap having the size de- 
pendence of AE( J2, h s , L) cx 1/L, the obtained value 
can be regarded as the L-dependent transition point, say 
h*(J2, L + 1). Then, extrapolating them to the thermo- 
dynamic limit, we estimate h*{J2)- Alternatively, as ex- 
plained in Refs. |2^ and |2j| — to which we refer the inter- 
ested readers for a detailed explanation — at h s = 0, the 
dimer excitation is in the subspace of total spin = 0, 
wave number k = ir, space inversion P = +1, and spin 
reversal T = +1 for L = (mod 4) [S% = 0, k = 0, 
P = -1, and T = —1 for L = 2 (mod 4)]. Thus, we can 
calculate AE(J 2l h s ,E) = E^(J 2 ,h s ,L) - E g (J 2 ,h s ,L) 
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from En(J2,h a , L) connecting to the dimer excitation and 
E s (J2,h s ,L), i.e., the ground-state energy. 
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FIG. 2. The boundary line of the Neel and "dimer" phases 
(for an explanation of "dimer" phase, see Ref. |l|). In- 
set shows some of extrapolations of the L-dependent criti- 
cal fields to the thermodynamic limit, where fitting curves 
are given. The double circle shows the fluid-dimer transi- 
tion point (J2,h s ) = (J|,0) at which the criticality of the 
boundary changes from the Gaussian to the Ising type. The 
spin-fluid (SF) state exists on the x axis of J2 < J|- 
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(a) Gaussian (v=2/3) 




(b), (c) Ising (v=l) 
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FIG. 3. The FSS plots of the excitation gap for systems 
of L — 20 — 28 in log- log scale. We have used v — 2/3 for 
the Gaussian transition [(a) J2 = JJ], and v — 1 for the 
Ising transition [(b) J2 =0.46, (c) J2 =0.5]. The slope of the 
dotted line is 1 showing the expected asymptotic behavior of 
the scaling functions. 

In Fig. [I], we demonstrate L and h s dependences of 
LAE( J2, h s , L) at several values of J2. Systems up to 



L = 28 sites are treated, where the Lanczos algorithm is 
used to obtain eigenvalues of the Hamiltonian in specified 
subspaces. We can see that the L dependence of the 
crossing point is almost absent for large value of J2 (near 
J 2 = 0.5), while it is visible for the small value case. 

After evaluating h*(J2, L + 1), we extrapolate them to 
L — > 00; here we assume the following formula: 

h* s (J 2 ,L) = h*{J 2 )+aL- 2 (l + bL-") (w>0), (10) 

where the b term stands for a correction to the leading 
one, and four parameters /i*(J2), a, 6, and u> are deter- 
mined according to the least-square-fitting condition. We 
used the data of L = 18 — 28, and extrapolated them to 
h*(J2) as shown in the inset of Fig. ^, where from bot- 
tom to top a series of data with fitting curves are given 
in the increasing order of J2. Consequently, we obtain 
the critical line h*(J2) as shown in the figure. The RG 
eigenvalue of the scaling variable is "almost zero" (i.e., 
marginal), while that of the staggered magnetic field is 
3/2 on the point (J2, h s ) = (J|, 0). The phase boundary 
is thus expected to behave as h*(J2) oc (J 2 — -TJ) 1 ' 5 ^; 
which agrees with a weak J2 dependence of the line near 
the point. 

At this stage, we shall check the criticalities by the up. 
of the FSS analysis in the vicinity of the boundary lines,EII 
i.e., we assume the following one-parameter scaling form 
of the excitation gaps in the finite system of L: 



AE{J 2) h s ,L) = L- l ^{L{h s - K{J 2 )) V ), 



(11) 



where v = 2/?> for J2 < J| and v = 1 for J2 > J% ■ Fur- 
ther, the asymptotic behaviors of the scaling function are 
expected as ty(x) oc x for large x and ^>(x) ~ const for 
x — > 0. Using the obtained transition points, we plot Eq. 
( pd| ) for the Gaussian (J2 = J^) and the Ising transitions 
(jf =0.46 and 0.5) in Fig. [| The results show that the 
data of different system sizes are collapsed on the sin- 
gle curve in both cases, and the asymptotic behaviors of 
>F agree with the expected ones (dotted lines) despite of 
the small L. From these plots, we can check the crossover 
behavior of the transitions driven by the staggered mag- 
netic field in the present frustrated quantum spin chain 
system and also the accuracy of the phase boundary line 
in Fig. H This good FSS behaviors may rely on the con- 
ditions that the marginal operator that brings about the 
multiplicative logarithmic correction to the pure power- 
law singularity of £ is absent on the fluid-dimer transition 
pointJHj and that the crossover region for J2 > 0.4 may 
be large enough to be detected. On the other hand, the 
FSS nature becomes obscure in the weak and intermedi- 
ate values of J2 since our system sizes are too small to 
reach the region. 

Finally, we evaluate the values of the central charge 
through the L dependence of the ground-state energyO 



E s (J 2 ,K(J 2 ),L)~e L 



7TVC 

~6L 



(12) 



(eo is the energy density in the thermodynamic limit). 
For this calculation, we should estimate v in advance; 
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here we use the FSS relation AE^ ~ (= \)/L and 

the least-square-fitting procedure for the data. The ob- 
tained results agree well with the Ising one for the con- 
siderably large values of J2, i.e., (v, c) ~ (1.439,0.494) at 
J 2 =0.46 and (1.248,0.491) at J 2 =0.50. For comparison, 
we also estimate them at J2 = 3 2 and obtain (v, c) ~ 
(2.365,0.995) [here the relation A£i ~ 2nvx 1 (= \)/L 
has been used to estimate v] . Consequently, we can again 
confirm the universalities of critical systems on the phase 
boundary line. On the other hand, we could not extract 
the reliable data using the above procedure in the small 
and intermediate regions of J2; this may be due to the 
finite-size effects on the line so that more detailed anal- 
ysis might be required in this region. 

To summarize, we have numerically investigated the 
ground-state phase diagram of the one-dimensional J\- 
Ji model in the staggered magnetic field. The crossover 
behavior of the second-order phase transitions driven by 
the staggered magnetic field occurs between the Gaus- 
sian and the Ising fixed points. According to the opera- 
tor correspondence between these fixed points and using 
the level-spectroscopy technique, we have analyzed the 
^2-odd excitation gap by the use of the phenomenolog- 
ical renormalization-group method, and determined the 
boundary line h*(J 2 ) representing the massless flow con- 
necting the fixed points. To check the two criticalities, 
we have performed the finite-size-scaling analysis of the 
excitation gap for typical parameter values, and we have 
also evaluated the values of the central charge. 

The present investigation has based upon the recent 
development of the (l+l)-dimcnsional double-frequency 
sine-Gordon theory, where its very interesting applica- 
tions in wide areas of researches including condensed 
matter physics have been pointed out.E3o However, as in 
the present case, numerical treatments of finite-size sys- 
tems may be required for quantitative discussions on the 
lattice Hamiltonian models. We think that our numeri- 
cal approach can also serve for the investigations of uio. 
complicated systems such as the correlated electrons £3i 
we will report on the application results of this approach 
in the future publications. 
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